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ABSTRACT 

We discuss one of the possible origins of large-scale magnetic fields based on a 
continuous distribution of toroidal electric current flowing in the inner region of the 
disc around a Kerr black hole (BH) in the framework of general relativity. It turns out 
that four types of configuration of the magnetic connection (MC) are generated, i.e., 
MC of the BH with the remote astrophysical load (MCHL), MC of the BH with the 
disc (MCHD), MC of the plunging region with the disc (MCPD) and MC of the inner 
and outer disc regions (MCDD). It turns out that the Blandford-Znajek (BZ) process 
can be regarded as one type of MC, i.e., MCHL. In addition, we propose a scenario 
for fitting the quasi-periodic oscillations in BH binaries based on MCDD associated 
with the magnetic reconnection. 
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1 INTRODUCTION 

According to observations, large-scale magnetic fields exist 
in many astronomical cases, such as galaxies, young stel- 
lar objects, neutron stars and black hole (BH) X-ray bina- 
ries (Han & Qiao 1994; Livio Ogilvie & Pringle 1999). It 
is widely believed that magnetic fields play an important 
role in dynamics of accretion discs and jet formation. How- 
ever, the origin of large-scale magnetic fields remains elusive 
(Blandford 2002). 

Recently, much attention has been paid to the magnetic 
connection (MC) of a rotating black hole (BH) with its sur- 
rounding accretion disc (Blandford 1999; Li 2000; Wang, 
Xiao & Lei 2002; Uzdensky 2004, 2005). Not long ago, Li 
(2002) calculated the magnetic field configuration of MC 
by assuming a single electric current flowing in equatorial 
plane of a BH. Wang et al. (2007) discussed the MC be- 
tween plunging region and disc (MCPD) as well as the MC 
between the BH horizon and the disc (MCHD). Very re- 
cently, Ge et al. (2008) analyzed the topology and feature of 
a magnetic field configuration arising from double counter 
oriented electric current-rings in an accretion disc around a 
Kerr BH. 

Enlightened by the above works, in this paper, we in- 
tend to investigate the magnetic field configuration arising 
from toroidal electric current distributed continuously in a 
thin disc around a Kerr BH, and discuss the MC of the 
BH with its surrounding disc. In addition, we propose a sce- 
nario for fitting the quasi-periodic oscillations in BH binaries 



based on the MC associated with magnetic reconnection in 
the disc. 

Throughout this paper the geometric units G = c = 1 
are used. 



2 MAGNETIC FIELD CONFIGURATION 
AND MAGNETIC FLUX 

In this paper, the magnetic field configuration is produced 
by the toroidal electric current flowing in the inner region of 
an accretion disc based on the following assumptions. 

(1) The disc is thin and Keplerian around a Kerr BH, 
and the concerned Kerr metric parameters are given by 
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(2) The toroidal electric current lies in the equatorial 
plane, varying continuously with the disc radius in a power- 
law, 



j = jo (r/rms) " 5{c 



^ r < Xr„ 
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where n is the power-law index for the variation of the 
electric current, rma is the innermost stable circular orbit 
(ISCO), and A is a parameter to adjust the radial width of 
the current distribution. 
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Thus the four-current vector flowing on the circle of 
radius r' is given by 



dip ' ''■='■ 



(3) 



Based on the above distribution of the toroidal elec- 
tric current we have the magnetic flux through a surface 
bounded by a circle with r = const and 9 = const as follows, 



* (a«, n, A, r, 6^,) = 27r ■ / dA^ (a,, n, A, r, 61; r') 



(4) 



In equation (|4]) dA^p {a, ,n, X,r,9;r') is the toroidal com- 
ponent of the electric vector potential, which can be deter- 
mined by the electric current loop located at r' — r' + dr' 
(Znajek 1978; Linet 1979), and the position {r,9) of the cir- 
cle is described in the spherical coordinates. 

A mapping relation between two circles with spherical 
coordinates {ri,9i) and {r2,92) is given based on the con- 
servation of magnetic flux as follows. 



*(a*,n, A,fi,6'i) = ^ {a,,n, X,f2,92) ■ 
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In equation (O is the dimensionless magnetic flux, 
being defined as 



* (a* , n. A, f , 9) = 2tvA^ {a*,n,X,f,9) / BqM'^ 
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where we have f = r/M and Bo = "2 jo- Based on (|6} we 
have the relation between the radii r^s and ro as follows. 



* (a*, n. A, fo,7r/2) = * (a*, n. A, fms,n/2) . 



(7) 



Incorporating equations ([T|) — ((Tj), we have the poloidal 
magnetic field configuration generated by continuous 
toroidal electric current with different values of the param- 
eters a«, n and A as shown in Fig.l. 

As shown in Fig.l, the four types of MC configurations 
can coexist for different values of a* , n and A, where the char- 
acteristic field lines of Types I, II, III and IV are indicated 
by lines 1, 2, 3 and 4, respectively. These MC configurations 
are described as follows. 

Type I: MC of the BH with the remote astrophysical 
load (MCHL); 

Type 11: MC of the BH with the disc (MCHD); 

Type III: MC of the plunging region with the disc 
(MCPD); 

Type IV: MC of the inner and outer regions of the disc 
(MCDD). 

Among the above configurations MCHD and MCPD 
have been discussed by a number of authors (Blandford 
1999; Krolik 1999; Gammie 1999; Li 2002; Wang et al. 2002; 
Wang et al. 2007), and these two MC processes are re- 
garded as the variants of the Blandford-Znajek (BZ) process 
(Blandford & Znajek 1977). On the other hand, the BZ pro- 
cess can be also regarded as one kind of MC, i.e., MCHL, 
which was formulated with the MC process in a united model 
(Wang, Xiao & Lei 2002). 

Based on equation ([7]) we have the magnetic flux 



)=ir/2 



penetrating the equatorial plane of a Kerr BH 



versus the disc radius f for different power-law index n as 
shown in Fig. 2. 

From Fig. 2 we obtain the following results. 

(1) The magnetic flux tjj{r,9)^^_ varies with f non- 
monotonically, and attains its peak value outside ISCO. The 
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Figure 1. The magnetic field configuration generated by toroidal 
electric current distributed continuously over the inner region of 
a thin disc around a Kerr BH with (a) a* = 0.5, n = 3.0, A = 3, 
(b) a, = 0.9, n = 3.0, A = 3, (c) a, = 0.5, n = 5.0, A = 3 and (d) 
a-t = 0.5, ra = 3.0, A = 5. 
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Figure 2. Magnetic flux penetrating tlie equatorial plane of a 
Kerr BH versus r for a, = 0.5 (a) A = 3 with n = 3,4 and 10, 
and (b) n = 4 with A = 2, 3 and 4. 



greater the power-law index n is, the less the peak value is, 
and the closer the position of the peak value is to ISCO. 

(2) Equation ((7|) provides a relation between the radii 
Vms and ro- For the given a*, the greater n corresponds to 
the less ro, and the greater A corresponds to the greater ro 
as shown in Figures 2a and 2b, respectively. 



3 FITTING QPOS OF BH BINARIES 

As is well known, kHz Quasi-Periodic Oscillations (QPOs) 
have been observed in several X-ray Binaries. As argued by 
van der Klis (2000), kHz QPOs in X-ray binaries probably 
originate from the inner edge of an accretion disc with a BH 
of stellar-mass, since millisecond is the natural timescale for 
accretion process in these regions. Although a number of 
models have been proposed to explain the QPOs in X-ray 
Binaries, no consensus has been reached on their physical 
origin (see a review by McClintock & Remillard 2006, here- 
after MR06). 



3.1 Fitting QPO frequencies of BH X-ray 
binaries based on MCDD 

Now we propose a scenario for fitting QPOs in BH binaries 
based on MCDD corresponding to the closed field lines indi- 
cated by line 4 in Fig.l. The field lines of MCDD are frozen 
in the disc, and the critical closed field line of MCDD con- 
nects the footpoint at ISCO and that at radius ro as shown 
by the thick solid line in Fig.l. 

The Keplerian angular velocity depends on the BH mass 
and spin, and decreases monotonically with the increasing 
disc radius, and it reads (Wang et al. 2003) 



(8) 

where uo is defined as i/q = (msif)~^3.23 X lO^Hz with 
tubh = M / Mq, and parameter ^ = r/rma is the disc radius 
in terms of r^s, and Xms = \frmIjM is a function of the 
BH spin a. (Novikov & Thorne 1973). 

Since the inner footpoints of the magnetic field lines ro- 
tate faster than the outer ones, the rotating disc will twist 
the field lines of MCDD, resulting in an increasing toroidal 
component of the magnetic field. This configuration allows 
antiparallel segments of the field lines to be brought in con- 
tact, giving rise to magnetic reconnection in a similar way 
to that between a central star and the accretion disc (Mont- 
merle et al. 2000). The magnetic energy thus librated could 
heat the plasma and ignites a fiare. In this way, successive 
reconnections and fiaring may continue periodically for some 
time, providing a possible way for producing QPOs in BH 
accretion disc. 

Considering that the magnetic reconnection arises from 
the twist of the closed field lines of MCDD, we infer that 
the flare occurs most probably in the critical field line, and 
the difference between the Keplerian frequencies at radii Vms 
and ro as follows. 



vmcdd 



nK(r„„)-nj^(ro) 



Xms +a*) 
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Thus we can fit the frequency vqpo of single- 
component QPOs as vmcdd based on equation ((9|, from 
which we find that the frequency vqpo depends on three 
parameters, ttibh, «» and ^o = ro/rms- By using equation 
((9]) we have an isosurface of vqpo— const in a parameter 
space consisting of these parameters. Taking the BH bi- 
nary GRO J1655-40 as an example, we have an isosurface 
of Vqpo = 300Hz with A = 3 as shown in Fig. 3. The pa- 
rameter 5o approaches the maximum 1.554 at point A with 
mpH ~ 6.6 and a, = 0.65, and the minimum 1.338 at point 
B with mpH ~ 6.0 and a» = 0.8. Thus we have n = 5.284 
and 6.656 by substituting ^0= 1.554 and 1.338 respectively 
into 0. 

We also have an isosurface of vqpo = 300Hz with a 
fixed value of n, e.g., n = 3, and thus we can determine 
the minimum and maximum values of ^o, and obtain the 
corresponding values of the parameter n. In this way we can 
fit the QPO frequencies of several sources based on MCDD 
as given in Table. 1. 



3.2 Improving the fittings for 3:2 QPO pairs 
based on MCDD 

As is well known, the 3:2 QPO pairs have been observed in 
the BH binaries listed in Table 1, i.e., GRO J1655-40, XTE 
J1550-564 and GRS 1915+105 (MR06). Also, the 3:2 QPO 
pair appears in the BH candidate H1743-322, and its be- 
havior resembles the BH binaries XTE J 1550-564 and GRO 
J1655-40 in many ways (Homan et al. 2005; Kalemci et al. 
2006; Remillard et al. 2002, 2006). 

Abramowicz & Kluzniak (2001) explained the pairs as a 
resonance between orbital and epicyclic motion of accreting 
matter. Recently, the resonance model is presented in a more 
realistic context, in which "parametric resonance" concept 
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Table 1. Fitting QPO frequencies of BH X-ray binaries based on MCDD. 





Input 






Output 










rriBH 


a. 


^max ) ^min 


A = 3 


n = 3 
X 




300 


6.6 


0.65 


1.554 


5.284 


1.547 


GRO J1655-40 


6.0 


0.8 


1.338 


6.656 


1.343 


450 


6.6 


0.65 


2.287 


3.551 


2.366 




6.0 


0.8 


1.628 


4.726 


1.649 




184 


10.8 


0.71 


1.485 


5.587 


1.483 


XTEJ1550-564 


8.4 


0.87 


1.233 


8.146 


1.242 


276 


10.8 


0.71 


2.036 


3.868 


2.083 




8.4 


0.87 


1.397 


5.778 


1.418 




113 


18.0 


0.98 


1.221 


7.122 


1.256 


GRS 1915+105 


10.0 


0.998 


1.114 


12.044 


1.119 


168 


18.0 


0.98 


1.360 


5.294 


1.424 




10.0 


0.998 


1.177 


8.714 


1.182 



Note. — The values of the BH mass and spin of the sources are adopted from RemilUard & McClintock (2006), Davis et al. (2006) and 
McClintock et al. (2006). The quantities ^max and ^min are the maximum and minimum of on the isosurface of vqpq = const. 




Figure 3. An isosurface of vqpq = 300Hz for GRO J1655-40 in 
the parameter space consisting of the parameters mgH,a, and 
Co- 



ls introduced to describe the oscillations rooted in fluid flow 
where there is a coupling between the radial and polar coor- 
dinate frequencies (Abramowicz et al. 2003; Kluzniak et al. 
2004; Torok et al. 2005). 

Not long ago, Aschenbach (2004, hereafter A04) found 
that the BH mass and spin are strongly constrained by the 
3:2 QPO pairs, and the ratio of the vertical epicyclic fre- 
quency to the radial epicyclic frequency is given by 

\ t7v(a.,f=32)/f^flK,r32) =3/2. ^ 

In equation (|10p Q,v and SIh are the vertical and ra- 
dial epicyclic frequencies, respectively (Nowak & Lehr 1998, 
Merloni et al. 1999). It has been found in A04 that equa- 
tion (|10p has only one solution for the two commensurable 
orbits, i.e., fai — 1.546 and r32 = 3.919 with an extremely 
high spin a, = 0.99616. As argued in A04, the radial gra- 
dient of the orbital velocity of a test particle changes sign 
in a narrow annular region when a* > 0.9953, and the BH 
mass of these binaries can be constrained in a very narrow 



range, which are consistent with the dynamically determined 
masses within their measurement uncertainty range. 

It has been pointed out that QPOs are generally asso- 
ciated with the steep power-law (SPL) state in BH X-ray 
binaries. Although the 3:2 QPO pairs could be interpreted 
in some epicyclic resonance models, there remain serious un- 
certainties as to whether epicyclic resonance could overcome 
the severe damping forces and emit X-rays with sufficient 
amplitude and coherence to produce the QPOs (e.g., see a 
review in MR06). The mechanism of producing QPOs based 
on MCDD might remedy the disadvantage of the epicyclic 
resonance models. 

The 3:2 QPO pair has been observed in H1743-322 
(Homan et al. 2005; Remillard et al. 2006; Kalemci et al. 
2006), and its behavior resembles the BH binaries XTE 
J1550-564 and GRO J1655-40 in many ways (Remillard et 
al. 2002, 2006). 

Based on the radial position of the vertical and ra- 
dial epicyclic resonance given by A04, we set fo — f32 
by adjusting the parameters n and A in equation ((TJ with 
a* = 0.99616. Thus we expect that both the vertical and ra- 
dial epicyclic resonance are stimulated and maintained due 
to energy transferred magnetically from the inner footpoint 
to the outer footpoint at fo. According to the strict con- 
straint required by the 3:2 QPO pair, we can estimate the 
BH mass as well as vmcdd as listed in Table. 2. 

In the above fitting the parameters n and A are not 
independent, being related by equation ((Tj for a, = 0.99616 
and fo = f32 = 3.919. And we have the curve of n versus 
A given in Fig. 4. Inspecting Fig. 4, we find that the power- 
law index n increases very sharply with the parameter A, as 
the latter is small, while n increases very slowly with A as 
the latter is great. These results imply that the fittings of 
QPOs based on MCDD are very sensitive to the parameter 
n, if the radial width of current distribution is small, while 
the fittings are almost independent of n, provided that the 
current distribution is wide enough. 

It is interesting to note that vmcdd for each source in 
Table 2. is about three times the corresponding upper fre- 
quency of the 3:2 QPO pair, i.e., vmcdd ~ 3i^up. Based 
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Table 2. Fitting QPOs in BH binaries based on the epicyclic resonance model. 



Source i^up/i^domrL rnsH a,, go vmcdd i^mcdd/'^up 



GRO J1655-40 450 / 300 



XTEJ1550-564 276 / 184 



GRS 1915+105 168 / 113 



H 1743-322 240 / 160 



6.77 1371.51 3.0478 



6.75 



11.06 



11.02 

■^j^g 0.99616, 2.9997 



17.77 



12.70 

12.66 733.421 3.0559 



1375.57 


3.0568 


839.521 


3.0417 


842.568 


3.0528 


502.169 


2.9891 


522.516 


3.1102 


731.111 


3.0463 



Note. - The BH masses of GRO J1655-40, XTEJ1550-564 and GRS 1915+105 are given by A04, and the the BH mass of H 1743-322 is 
determined by the same way given in A04. The parameter go = 2.9997 is determined by go = = To/rms for a* = 0.99616. 



3.5 
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A 

Figure 4. The curve of the parameters n versus A constrained 
by the 3:2 QPO pairs of the BH binaries. 



on general relativity we have vmcdd / Vup = 3.05, which is 
independent of the BH mass. This result seems consistent 
with the commensurate frequencies of the 3:2 QPO pair. In 
addition, the combination of MCDD with the epicyclic res- 
onance might be helpful to interpret the association of the 
QPOs with the SPL states in BH binaries, because both 
the steep power-law component of radiation and the QPOs 
in BH binaries have the same origin of the magnetic recon- 
nection. We shall address this issue in details in our future 
work. 
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